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PACS 03 . 75 . Gg - Entanglement and decoherence in Bose-Einstein condensates 

Abstract - The evolution of an interacting two-component Bose-Einstein condensate from an 
initial phase state leads to a spin squeezed state that may be used in atomic clocks to increase the 
signal-to-noise ratio, opening the way to quantum metrology. The efficiency of spin squeezing is 
limited by the finite temperature of the gas, as was shown theoretically in a spatially homogeneous 
system. Here we determine the limit of spin squeezing in the realistic trapped case, with classical 
field simulations, and with a completely analytical treatment that includes the quantum case. 
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Introduction. — In good atomic clocks, the signal- 
to-noise ratio is determined by the quantum noise, that is 
^the partition noise of uncorrelated atoms among two in- 
Olfernal levels a and b, rather than by technical noise Q~|. 

'i'he resulting statistical uncertainty on the transition fre- 
^HQuency w a & is then Acj Q b = l/(riV 1 / 2 ) where N is the 

atom number and r the Ramsey interrogation time. 
G\ ' One can beat this so-called "standard quantum limit" 
by introducing appropriate quantum correlations among 
i^^the atoms, that is using squeezed states of the collective 
•spin [2] having reduced fluctuations AS± along some di- 
^^rjection transverse to the mean spin (S). The resulting 
statistically uncertainty on the transition frequency is re- 
i— I duced by the squeezing parameter £ < 1 to the value [3] 
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A fundamental issue is then to know how £ scales in 
the large- N limit. Practical realisation of spin squeezing 
requires some non-linearity in the spin dynamics [H[S|. In 
two-component Bose-Einstein condensates, it is provided 
by the atomic s-wave interactions [6], and first experi- 
mental implementations have been performed with a dy- 
namical control of the interactions [7] . The original model 
Hamiltonian proposed in [5] is then realized if one restricts 
to the two condensate spatial modes (in a and in 6), re- 
sulting in a squeezing parameter minimised over time such 
that £ 2 lin ~ N- 2 ? 3 . For one million atoms, this predicts 
a signal-to-noise increase by a factor 100. There may be 
however serious shortcomings of the two- mode model. The 
combined effect of one-body and three-body particle losses 
was shown in [5] to impose a non-zero lower bound on £ m j n . 



More fundamentally, the finite temperature gases in the 
experiments are intrinsically multimode; in the spatially 
homogeneous case, this was predicted in [9j[T0] to lead to 
a non-zero value of in the thermodynamic limit, in 
sharp contrast with the N~ 2 / 3 scaling, as a consequence 
of random dephasing of the condensate introduced by the 
non-condensed quasi-particles [llj . 

Here, we perform the last step towards a full theoretical 
understanding of the spin squeezing dynamics, by includ- 
ing the effect of the trap present in the experiments. Al- 
though it significantly complicates the theoretical analysis, 
we will see that an analytical treatment can be pushed to 
the end to obtain £ m i n in the large system-size limit. 

Classical fields. — Classical field models require an 
energy cut-off of the order of fcgT, where T is the initial 
temperature of the gas. Here it is provided by a cubic 
lattice model, with a unit cell volume dV adjusted to re- 
produce the quantum non-condensed particle number for a 
zero chemical potential spatially homogeneous ideal Bose 
gas. The corresponding Hamiltonian is 



H = E T, dV 



a—a.b 
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with a one-body part ho = p 2 /(2m) + U(r) involving 
the kinetic energy part [eigenmodes are plane waves of 
wavevector k in the first Brillouin zone] and the internal- 
state independent trapping part U (r) . We assume that 
the coupling constant g is the same in the two internal 
states a and b, as for F = 1 spinor condensates, and that 
there is no interaction between a and b, which may be 
realized by spatial separation of the two components. Ini- 
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tially, at time t = 0~, the gas is in thermal equilibrium 
in the internal state a, with a fixed total particle num- 
ber iV and a coupling constant g(0~). To mimic quantum 
noise in the initially empty internal state 6, we sample 
the Gaussian Wigncr distribution of the vacuum state in 
b, which amounts to introducing on average 1/2 fictitious 
particles in each of the M spatial modetlj- At t = 0, the 
gas experiences a tt/2 pulse, inducing the field mapping 



V>a,&(r, 0+) 
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This distributes on the average half of the particles in each 
internal state, with Poissonian fluctuations, and it puts the 
system out of equilibrium. To minimize the excitation of 
a breathing mode of the condensate [14] , we subsequently 
increase the coupling constant by a factor two, 



(4) 



ff (0+) = g = 2g(0-), 



which may be realized using a Feshbach resonance □ Q 
The before-pulse coupling constant <?(0~) is related to 
the before-pulse s-wave scattering length by <?(CP) = 
Airh 2 a(Q~)/m. After the pulse, each field ip a ,b is evolved 
according to the non-linear Schrodingcr equation 



ihdtipo = [h + g\i'a(r,t)\ 2 ]tp a 



(5) 



The variances of the various collective spin components 
S x + iS y = Er dV ^*a( r )Mr) and S z = (N a - N b )/2, 
where N a is the particle number in internal state a, give 
access to the spin squeezing parameter £ 2 (i) as a function 
of time. According to Eq. (UJ), one needs in particular [§]: 
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where (. . .) represents the average over all stochastic re- 
alisations. In practice, spin squeezing is supposed to be 
applied to large atomic ensembles, as in clocks, so we con- 
centrate our study on the thermodynamic limit, where the 
single-particle level spacing tends to zero for fixed temper- 
ature and chemical potential. We then find as in [51 [TO] 
that £ 2 (i) is a very flat function around its minimum, and 
that its minimum tends to a finite value in the ther- 
modynamic limit, see inset in Fig. \T\ This is a striking 
consequence of the multimode nature of the problem. 

In the spatially homogeneous case, in the weakly inter- 
acting regime, £^j n obeys a simple scaling formula [51 ITU]. 
Its intuitive generalisation to the trapped case is 
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1 The Wigner method in the context of spin squeezing was also 
used in 1121 and recently extended to the case with particle losses 

□i. 

2 Strictly speaking, one should take 9(0")^ = g{0+)(N +Af)/2. 
Here, however, one has N 3> Af/2 as required in 1151 . 

3 We have however seen numerically that for g(0+) = g(0~) such 
excitation does not have a large impact on the squeezing. 
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Fig. 1: Minimal spin squeezing parameter £ min in an isotropic 
harmonic trap, U(r) = mu) 2 r 2 /2, as a function of temperature. 
After rescaling as in Eq. 0, it is universal in the thermody- 
namic and weakly interacting limit. Analytical result (|33|l for 
quantum fields: solid line. Classical field model: simulations 
(squares) vs analytics (|26p (crosses). Inset: £ 2 (t) from simu- 
lations approaching the thermodynamic limit; feT/(i = 2.89, 
[p(0,0")a 3 (0~)] 1/2 = 1.4 x 10" 4 , N from 10 5 to 5.6 x 10 5 . 



where fi is the Gross-Pitaevskii chemical potential of the 
gas, and p(0,0~) is the mean atomic density before the 
pulse in the trap center. We have checked numerically, 
by producing data with various atom numbers and inter- 
action strength, that this scaling holds for the trapped 
system. The rescaled minimal spin squeezing £^ in as a 
function of the rescaled temperature, that is the function 
/, is plotted with symbols in Fig. [TJ The temperature is 
limited to the range ksT > /i where our classical field 
model makes sense. 

Analytical approach for classical fields. We 

develop analytics within the framework of the number- 
conserving Bogoliubov approach [TrjlITT] , along the lines of 
our studies of the homogeneous case [TU] . The Bogoliubov 
approach is used first to describe the initial distribution 
of the fields (before the pulse); this is relatively standard, 
and we shall be brief on this aspect here. Second, it is used 
to describe the dynamical evolution of the fields after the 
pulse, which is more involved as we now explain. 

The theory introduces as a small parameter the non- 
condensed fraction of the gas, £Bog, and it is pushed 
to first order included in eB 0g - It would remain quite 
involved without the use of a second small parameter, 
£sizo = 1/A 1 / 2 where N is the total particle number, con- 
trolling the approach of the thermodynamic limit Q In 
the Bogoliubov approach, one singles out the contribution 
of the macroscopically populated modes, corresponding to 
the condensate wavefunctions 4> a in each internal state: 



Vv(r,i) = a a (t)(f> a (r,t) +tp± a (r,t) 



(8) 



4 The two small parameters are not fully independent. Since the 
system size is 2> the thermal wavelength, one finds Nesog 2> 1- 
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The non-condensed fields ip±a are orthogonal to the con- 
densate wavefunctions; their squared modulus gives the 

non-condensed density so that they are smaller than the 

1 li 

condensate fields a CT </> CT by a factor e Bog - This allows the 
Bogoliubov expansion of the Hamiltonian in powers of ip± . 

For each component a the condensate wavefunction </v 
is an averaged quantity: it is the most populated cigen- 
state of the one-body density matrix (ip*(r',t)ip a (r,t)}. 
Due to symmetry in our problem, one has cj> a = </V The 

1/2 

expansion of CT in powers of e B reads, with notations 
slightly different from [T7] : 



(9) 



The zeroth order contribution $ corresponds to Gross- 
Pitaevskii theory; it is time-independent here, due to the 
trick ((4]), so it can be taken as the real and positive solution 
of the stationary Gross-Pitaevskii equation 



0= [h + giNv}® 2 ~ ^ 



(10) 



where (N a ) is the mean particle number in each internal 
state after the pulse. The equation for the correction (jP^ , 
due to condensate depletion and interaction with the non- 
condensed particles, is not needed; wc shall use that 4>^ 
scales as £Bo g ^, and that it is orthogonal to 3> [T7] . 

Of crucial importance is the modulus-phase representa- 
tion of the condensate amplitudes 



a<r(t) 
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where the condensate phases a and particle numbers n a 
are fluctuating and timc-dcpcndcnt quantities. First, in- 
clusion of this phase factor in the non-condensed fields 
leads to the fields (called A ox in [T7]) that are the main 
dynamical variables of the number-conserving theories: 

A ff (rJ) = e ^W^(r,t) (12) 

Second, the phases 9 a {t) are unbounded variables that in- 
crease linearly in time, contrarily to the condensate num- 
bers n a (t) or the fields A^- that merely oscillate in time. 
The simulations show that the best squeezing time remains 
finite in the thermodynamic limit, as it was the case for 
the spatially homogeneous system [MTU] ; the phase differ- 
ence {0 a — 9 b )(t) thus remains Oil/N 1 / 2 ) and <C 1 over the 
relevant time window. Linearizing the quantity e l ^ 0b ~ ea ^ 
that appears in the expression of the collective spin com- 
ponent S y , as in [9ll0j , one obtains a contribution growing 
linearly in time that eventually dominates over the collec- 
tive spin component S z (which is a constant of motion). 
The resulting squeezing parameter £ 2 (i) reaches its mini- 
mal value at large times |£o 
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{{e a -6 b ){N a -N b )Y 



large t ((0 a - 6 b ) 2 )({N a - N b )*) large t 



Cin (13) 



5 As in the homogeneous case, large times correspond (in the 
Thomas-Fermi regime) to £ m i n /^f/?l 2> 1. 

6 As in 1101 we use the fact that, at large times, S y ~ {Of, — Oa)^ 
where T = ^/n a nf, + Re z >^ r dVA*A b can be replaced by its mean 
value (J 7 ) in the thermodynamic limit since 9 a — b is of zero mean. 



Note that N a — N b scales as iV 1 / 2 , since it is of zero mean 
and of variance ~ N. 

We thus need to evaluate the time evolution of the con- 
densate phases. From the equation of motion ([5]), one gets 
the expression of 4jQ<t that we transform with the splitting 
([5]). keeping terms up to order one included in eBog and 
up to order one included in e s i zo , and using where neces- 
sary n<7 = N a — '^2 r dV\K a \ 2 . Going beyond first order 
in eBog is indeed beyond Bogoliubov theory. According to 
(|13|) . the relevant quantity is N 1 ^ 2 (6 a — 6b), and terms in 
the phase difference that are o(e s i ze ) are negligible in the 
thermodynamic limit. Also, the crossed terms scaling as 
r e S j ZC are negligible. We finally obtain after the pulse: 



- h-9 a ~ n + (N a )g J $ 3 (0^ + c.c.) 
+ (n a -(N a ))g f^ + 2g /<I> 2 (W + 



+ ((^) 1/2 + fc)j * 3 (A„ + A:), (14) 



where J is a short-hand notation for ^ r dV and we have 
introduced the opposite quantities of order unity cor- 
responding to the fluctuations in the particle numbers, 
q a = {N a - (N a ))/(N a y/ 2 . The first line in (0 is the 
same for the two internal states due to symmetry and will 
cancel out in the phase difference 8 a — 8 b . The second line 
was already present in the spatially homogeneous case [TS] . 
The last term in (|14p is new; it is non-zero, even if A CT is 
orthogonal to </> CT , because <& and <f> 3 are not proportional 
in a trap. This new term is linear in the non-condensed 

1 /2 

field and thus potentially larger by a factor l/e Bog than 

the preceding term. It forces us to evaluate the fields A^ 

1/2 

one order in e Bog beyond the usual Bogoliubov approxima- 
tion, which makes the trapped case treatment more subtle 
than the homogeneous case one, as we now see. 

In the usual number-conserving Bogoliubov approaches, 

1 /2 

the fields A CT are evaluated to the leading e B ' og order only. 

They are expanded on the Bogoliubov mode functions 

Ufc(r) and i>fc(r), that are orthogonal to $ and normalized 

as (uk\uk) — (vk\i>k) = 1- They are here common to both 

internal states and most importantly, thanks to the trick 

l[4|). they are at t < and t > the same time-independent 

eigenstates of the Bogoliubov operators of eigenencrgics 

1/2 

ej; > 0. To first order included in e BoE , the modal expan- 
sion then reads (after the pulse) 

A«(r,t) =^cU0 + )e- iekt/h u k (r)+cU0+)e i ^%i(v) 
fc 

(15) 



From ([3]) and along the lines of [TD], one finds 
c fcCT (0+) = -j=[ Cka (p-)^B k ], 



(16) 



omitting a (^(l/iV 1 / 2 ) contribution involving the conden- 
sate phase change due the pulse, negligible in the large 
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N limit [TO]. The modal coefficients Cfc a (0~) in the ini- 
tial thermal equilibrium state have independent Gaus- 
sian distributions, with (|cfc a (0 _ )| 2 ) = fc^T/efc according 
to the equipartition formula. The amplitudes Bk result 
from the projection of the initial Wigner noise in internal 
state b over the Bogoliubov modes, and are thus Gaussian 
variables of zero mean statistically independent from the 



Cfe a (0 )'s. We only need here that (\Bk 



5 + (vk\vk)- 



Due to the potentially much larger last term in Eq. (fT4|) . 
the leading order (|f 51) is not accurate enough, and one 

1/2 

must go to the next order in e Bog - Furthermore the fluctu- 
ations of the particle numbers N a will give a contribution 
not included in (|15p . since the usual number-conserving 
theory assumes a fixed total particle number. One writes 
A CT = e~ l6 "Q < j > ip a , where Q§ projects orthogonally to <f) a , 
and one calculates the temporal derivative. After inser- 
tion of the splitting © and of the expansion © , one gets 

■ such terms contribute to the 



terms scaling as fiA a e% og e^ 

last term of Eq. ([T4f as M e Bog e size! so tncy navc to be 
neglected if a + 1/2 > 1 or /3 > 1 or if a + 1/2 = /} = 1. 
This leads to 

ihdtK Q(h Q — fi + 2(A (T ) 5 $ 2 )A ff + Q(N a )g$ 2 A* a 
+ (n a -{N a )) (N^Qg® 3 + q <J {N a ) 1 ' 2 [Qg<5> 2 (2 A a + A; ) 



gA„ * 4 -]+(N t ,)V 2 [{N a )Qg& 



1 / 2 3 g{$(2|A (T | 2 + a 2 ) - k( [$ 3 (a ct + a;)]} 



(17) 



where Q projects orthogonally to $. The various terms 
in the right-hand side have a simple physical interpreta- 
tion. The first two terms constitute the usual Bogoliubov 
equations of motion, written for a number (N^) of par- 
ticles, and whose are the positive-energy eigen- 
modes. Due to the fluctuations of the particle numbers, 
these usual Bogoliubov equations differ from the ones writ- 
ten for a given realisation with N a particles, hence the 
fourth term (on two lines) . The third term corresponds to 
a similar effect at the level of the Gross-Pitaevskii equa- 
tion, further including the condensate depletion: It takes 
into account the fact that the Gross-Pitaevskii solution for 
n a particles in the condensate of internal state a slightly 
differs from the solution $ for (N a ) particles in that con- 
densate, providing a source term for A a . This third term 
already appeared in Eq. (3.14) of [19] (except for the miss- 
ing projector Q). The fifth term, involving the deviation 
(f>( 2 > of the condensate wavefunction from gives the same 
contribution for a = a,b; to the relevant order, its con- 
tribution will cancel out in the phase difference 8 a — Of,. 
The sixth term describes the effect of the cubic interaction 
among the Bogoliubov quasi-particles; whereas this inter- 
action usually leads to a true transfer of quasi-particles 
among the Bogoliubov modes (according to the so-called 
Bcliaev-Landau processes) , its relevant effect here is rather 
a reactive effect that shifts the value of the field A„. 



We now apparently have to solve Eq. (fTT)) . This 
formidable task is greatly simplified by the fact that we 
need here, in the long time limit, only the linearly diverg- 
ing part of the phases 9 a {t). What should be evaluated is 
thus only the time averaged part Q a \ the oscillating parts 
9 a (t) — 9 a will be neglected in the so-called secular approx- 
imation for the condensate phase Q2] carefully justified for 
spin squeezing in the homogeneous case in [10] . In the sec- 
ond line of (fT4|) . the leading value (fT5|) for A CT suffices. In 
the third line of ((T4)) . it would give the insufficient result 
A CT = 0, so we have to determine the leading order non- 
zero approximation for A CT from (|17[) . In the subleading 
terms of ([XT]) , that is the fourth, sixth and seventh terms, 
we can replace A a by the leading value ([T5]l ; after tempo- 
ral average, the fourth term thus disappears. Adding the 
temporally averaged (|17|) to its complex conjugate gives 
the inhomogeneous equation 

= M(A a + A*) + S (18) 

with the hcrmitian operator M given by 

M = Q{h - /i + 3(A CT ).g$ 2 )Q (19) 

and an easy-to-reconstruct source S term that we shall not 
write explicitly. Since $ is the ground energy solution of 
the Gross-Pitaevskii equation, the operator M is positive 
and invertiblc (in the subspace orthogonal to $). Whereas 
the explicit expressions of M _1 , and thus of M~ 1 S, are 
unknown, in order to get the last line of (|14[) we fortunately 
only need 

M-^Qgtf = -d {Nr) *, (20) 

as can be deduced from the derivative of (fTU)) with respect 
to the particle number [T7]. Then we get 



ff $ 3 (A CT +A;)= (gQ<f> 3 )M-\-S) 



All this leads to 



(21) 



t- l m a {t) 



large t 



2|A 



(1)|2 



A 
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d {N ^((N a )g^ 2 ) 



- g{N a ) yy(Ai 1) (r)+c.c.)(Ai 1) (r') + c.c.)$ 3 (r)9 (Arff> <I>(r') 

+ terms symmetric under exchange a O b (22) 

where we used d^ N ^fj,— j g$ 4 = 2(N a )g j $ 3 9 (Ar<T )$. The 
last line in (|2"2l originates from the correction </>( 2 ) to the 
condensate wavefunction; as it assumes the same value in 
both internal states a and b, it cancels out in the phase 
difference, so that its explicit value is not needed here. 
Remarkably, the A-dependent terms in (f2"2"]l exactly re- 
produce the classical field version of the ones appearing 
in the expression of the full (not Gross-Pitaevskii) chem- 
ical potential /itot of the gas at thermal equilibrium, as 



p-4 



Limit of spin squeezing in trapped Bose-Einstein condensates 



given by Eq. (103) of [50]. This confirms the expectation 
that the condensate phase evolves at the average pulsation 
—Utot/fy in other words, this constitutes a test of (|22p. 

It remains to use (|13p to obtain the squeezing parameter 
minimized over time. As in our previous studies, the phase 
difference has the structure 



h{t)-e h {t) 



large t 



-[N a -N b + D]^t 



(23) 



where D collects the contributions that would be absent 
in the two- mode model realizing the proposal of |2]. The 
time average in (|2"2"|) suppresses the crossed terms that in- 
volve two different Bogoliubov modes, since they oscillate 
in time. This leads to a single sum over modes: 



D 



E*[i 



Cka\ 



C kb 



(24) 

where we used (JTBJ) to transform |cfe a (0 + )| 2 — |c£b(0 + )| 2 
and the fact that the quasi-particle occupation numbers 
|cfco_(t)| 2 are constants of motion within Bogoliubov the- 
orju. The general expressions for dk are easily deduced 
from (|22p . They look involved, but an inspired applica- 
tion of the Hcllmann-Fcynman theorem to the Bogoliubov 
operator [of which (uk,Vk) is an eigenmode] gives 



d {N a )£k 



d, 



N„)V 



(25) 



The simplicity of the result is understood from another 
expression of the gas full chemical potential /Lttot at ther- 
mal equilibrium, which is most rapidly obtained within 
the microcanonical ensemble as in |18j : /Lttot is then the 
derivative of the energy E to t with respect to the particle 
number N a at fixed entropy; since within Bogoliubov the- 
ory, E to t = E + J2k e k(\cka\ 2 ), and the occupation num- 
bers (|cfe CT | 2 ) are constant at fixed entropy, we get (|2~5"j) . 

Finally, inserting the structure (J23J) in f| 13[) and keep- 
ing the leading (that is first) order in eBog gives for the 
classical field modcH: 



£ 2 



class (D 



N 



= ^£4^(i + 2Kk>) 



(26) 



This expression is successfully compared to the numeri- 
cal simulations (see symbols in Fig. [T]), after a numerical 
diagonalisation of the Bogoliubov operator to obtain the 
cigenfunctions (uk,Vk) and the cigenencrgies ek- 

Extension to quantum fields. — The previous ana- 
lytical developments have been performed within the clas- 
sical field model, in order to allow for a comparison with 
the numerical simulations. They can be quite directly 



7 This requires that the squeezing time is shorter than the ther- 
malisation time. This holds in the weakly interacting limit [9], 

8 In particular, we used that {(N a —N),)D) ?s JVesog, see Appendix 
D in I10| , so that its contribution, which appears squared in the 
numerator of (1131) - is negligible as compared to (D 2 ) a NtBog- 



transposed to the quantum case, where the number con- 
serving Bogoliubov theories were initially developed. The 
condensate phases and particle numbers are now conju- 
gate hermitian operators [tv,^/] = i5 aiT i, and the classi- 
cal fields become operators ip a and A CT . Also Ck a ,Bk and 
D have quantum counterpart^. The energy cut-off is no 
longer restricted to ksT and one can take the limit of 
a vanishing lattice spacing in the Bogoliubov results for 
^in, so that the notation J reduces to a true spatial inte- 
gral and one at last studies the case of trapped particles 
in real continuous space. We obtain 



quant (Z) 2 



N 



k 



1 + 2<^ fc |^ fc ) 

e c k /k B T _ 1 



\Vk\Vk) 



(27) 

The coefficients dk have the same expression in terms of 
the Bogoliubov modes (uk, Vk) as in the classical field case, 
so that Eq. (|2"5|) still holds. As compared to (f26"P . the oc- 
cupation numbers of the Bogoliubov modes are now given 
by the Bose formula. There is also a quantum term that 
subsists at T = 0; the fact that £^ in is non-zero for T = 
is due to the excitation of quasi-particles by the 7r/2 pulse. 

Continuous spectrum limit. — Up to now, we have 
taken advantage of the large-iV limit (through the e s i zc 
expansion) but we have not explicitly used the fact that 
the thermodynamic limit also corresponds to a vanishing 
single-particle level spacing, where the spectral sum in the 
quantum result ([27]) may be replaced by an integral. This 
we now fully implement for an isotropic harmonic trap 
U(r) = muj 2 r 2 /2, where the semi-classical limit a la WKB 
corresponds to Huj/[1 — > and can reasonably easily be 
taken thanks to integrability of the classical motion issued 
from the Bogoliubov equations [3T] . 

For the condensate wavefunction, one takes the 
Thomas- Fermi approximation W(r) = (N a )g^ 2 (r) ~ 
hi - U(r)]Y\p - U{r)} [22], where Y is the Hcaviside 
function; $ then strictly vanishes at distances larger than 
the Thomas-Fermi radius i?, and from the usual Thomas- 
Fermi formula one obtains 



(fi/hw) 3 = 157V(V8) 1/2 [p(O,0-)a 3 (0-)] 1/2 . 



(28) 



Due to rotational symmetry, the Bogoliubov modes are 
labeled by the radial n > 0, angular momentum I > and 
azimuthal mi quantum numbers, so k = (n,l,mi). The 
mi -independent eigenenergies e n ,i, thus of 21 + 1 degener- 
acy, are approximated by the Bohr-Sommerfeld quantiza- 
tion rule [2T| 



Tih{n + 1/2) 



drp r (r, e fc ,/i) 



(29) 



9 In particular, one has [Bfc,BL] = Skk' an d (B^Bfc) = (yk\ v k)- 
10 Modcs of angular momentum I are then degenerate; their crossed 
terms in (122 b . not suppressed by time average, are killed by spatial 
integration, so that the form (1241) holds. 
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where the positive classical radial moment p r , such that 



p 2 r (r,E,n) = 2 + w21/2 
2m 



h 2 i 2 

2mr 2 



2W + U - fi 



(30) 

vanishes at the two classical turning points n and r 2 . Tak- 
ing the derivative of (|29|) with respect to [i gives 



drdftPrfaeicfj), (31) 



takes a finite optimal value in the thermodynamic and 
weakly interacting limit, a value that we determine ana- 
lytically. Our theory applies both for classical fields ([25)1 
and quantum fields f|2T[) . and for arbitrary trap geometries. 
It indicates that a large metrologic gain is still possible at 
finite temperature in realistic experimental conditions: for 
k B T ~ fj, ~ lO/to;, and assuming that one million is a large 
enough atom number to reach the thermodynamic limit, 
we predict a signal-to-noise increase by a factor l/£ ~ 30. 



where the angular frequency oj c \ of the classical mo- 
tion is given as usual in classical mechanics by tt/u> c \ = 
/J" 2 drdEPr(r,£k, [!■)■ The mode functions Ufc(r) and i/fc(r) 
have their angular part proportional to spherical harmon- 
ics Y™ 1 (9, 4>)', the WKB theory gives explicit expressions 
for the radial parts [H| leading to 



l + 2(vk\v k ) 



dr — 

n Pr 



(32) 



To evaluate the resulting integrals, which can be done 
analytically, it is convenient to introduce the reduced 
energy e = E//J, and the rescaled angular momentum 

j = inu/(2n). 

Two types of classical orbits are obtained in the trap 
[2"T] . A first type corresponds to orbits that are purely out 
of the condensate, due to a large enough angular momen- 
tum (and centrifugal barrier), and a not too high energy, 
1 < 2j— 1 < e < j 2 ; in this case, the quasi-particlcs simply 
experience the effective potential U — fx, so that dk — — 1, 
u c i = 2lj and (vk\vk) — 0. The second type corresponds 
to mixed orbits that cross the condensate boundary, that 
is ri < R < r 2 , which corresponds to the parameter range 
< j 2 < e. Also in that case dk is negative (at variance 
with the homogeneous case), but it is larger than — 1. This 
implies that £^ in is always smaller than the non-condensed 
fraction as in [9] . Integration over the purely external or- 
bits can be performed explicitly: 



2 

min 



[p(0,0-) O 3(0-)]Va 



i5y^ n B r 

2\/2 V M 



93 (e 



-fi/k B T 



)+/n 



(33) 

where g a (z) — J2 n >i zn / na ^ s ^ ne Bose function. The 
mixed-orbit contribution still involves a double integral, 
which is ultraviolet and infrared convergent: 



fn 



15V2 



de 



jdj 



J 



K , en 
— coth — — 
J 2k B T 



^ (34) 

where the functions /, J, K of e and j are given in the 
appcndbJ"! This leads to the quantum result plotted as 
a solid line in Fig. [TJ 

Conclusion. — We show that spin squeezing driven 
by interactions in trapped bimodal atomic condensates 



11 For kgT 3> fJ. external orbits dominate over /„ 
0(k B T/fi) 2 . 
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Appendix. — With the change of variable l—r 2 /R 2 = 



e sinh z, and then X 



3 ±2 



the semi-classical integrals can 



be evaluated. Setting J = -klj/uj c \ = J1 + J2, where J\ (J2) 
is the in(out)-condensate contribution, one finds: 



Ji = 
■h = 
Setting / = 



z/V2 



arccos 



2J 1 



1 1-e 
-arccos ^—^——^ 



1)V2 



= —Jd^ek = I\ + I2, one finds 1 2 = J2 and 

/l = 7= arccosh ——5 — — 

y/2 (2j 2 + e 2 + I) 1 / 2 



(35) 
(36) 

(37) 



Setting K = J(l + 2(v k \v k )) = K 1 +K 2 , one finds K 2 = J 2 
and 

* 1 = £( *t v/ !? ) -£ (38) 



2(e 2 + 2j2) 
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